Strong phonon-plasmon coupled modes in the graphene/silicon carbide heterosystem 
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We report on strong coupling of the charge carrier plasmon uj pl in graphene with the surface 
optical phonon ojso of the underlying SiC(OOOl) substrate with low electron concentration (n = 
1.2 x 10 15 cm -3 ) in the long wavelength limit (qn — > 0). Energy dependent energy-loss spectra 
give for the first time clear evidence of two coupled phonon-plasmon modes u>± separated by a gap 
between ujso - > 0) and uito (<7|| >> 0), the transverse optical phonon mode, with a Fano- 
type shape, in particular for higher primary electron energies (.Eo > 20eV). A simplified model 
based on dielectric theory is able to simulate our energy - loss spectra as well as the dispersion 
of the two coupled phonon-plasmon modes w±. In contrast, Liu and Willis [H] postulate in their 
recent publication no gap and a discontinuous dispersion curve with a one-peak structure from their 
energy-loss data. 
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The graphene silicon carbide heterosystem is a promis- 
ing system for the future application of graphene in 
micro- and nanoelectronics @, HI- Silicon carbide as a 
substrate for microelectronics is already used industrially 
and the epitaxial growth of graphene on silicon carbide 
has already been investigated for several years now 0, H| , 
and perfectionalized towards wafer scale homogeneous 
graphene [j| . Still, many of the interactions between the 
graphene and the silicon carbide substrate have yet to be 
understood. For example, the carrier dynamics may be 
strongly influenced by the long-range coupling to the po- 
lar modes of the substrate [fj( , which possibly results in 
a strong reduction of the graphene mobility, if compared 
to free standing graphene. This remote scattering can be 
important in future graphene devices. 

In this contribution we report about our experimen- 
tal investigation of the carriers in the conduction chan- 
nel with the long-range polarization field created at the 
conductor /dielectric interface. Emphasis is also given to 
the theoretical interpretation of the experimental inelas- 
tic electron scattering results by calculating the dielectric 
surface loss function. The coupling of collective electron 
(or hole) modes with optical phonons in semiconductors 
(e.g. InN, InP, GaAs and others) has already been a 
target of extensive investigations and helped_to under- 
stand important interface characteristics 0-13 . Unlike 
these conventional two dimensional electron gas systems 
(2DEG), graphene exhibits a linear electron dispersion 
relation, but the plasmon dispersion remains [10 . Fur- 
thermore, the almost vanishing damping of the plasmon 
mode and the strong spatial confinement, in contrast to 



other sheet plasmons observed so far [llj, [12| , makes it 
a showcase model for the investigation of the coupled 
phonon plasmon modes. 

Inelastic electron scattering utilizing special high res- 
olution monochromators and analyzers, also known 
as high resolution electron energy loss spectroscopy 



(HREELS), is used in surface science to investigate in- 
tentional and unintentional adsorbates as well as surface 
phonon and plasmon modes on a wide variety of materi- 
als [13(1 . Dispersion measurements can be obtained in dif- 
ferent measurement methods. By changing the analyzer 
rotation angle impact scattered electrons can be analyzed 
over the whole Brillouin zone, which is mostly employed 
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FIG. 1: Energy-loss spectra in specular direction at a spec- 
trometer resolution of 20 cm -1 on different samples, (a) 
on silicon carbide hydrogen etched with Eo = 36 eV, (b) 
on the graphene/SiC(0001) heterosystem with Eo = 20 eV. 
T = 300K. 
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for measuring the phonon- and electron-dispersion[14j . 
To investigate the dispersion very close to the center of 
the Brillouin zone the analyzer is kept in specular di- 
rection and only the primary energy of the impinging 
electrons is varied. 

We took all HREELS spectra in specular direction, 
with the impact and scattering angle 9 fixed at 64° rel- 
ative to the surface normal. The momentum transfer 
parallel to the surface for the dispersion measurements 
is calculated in this particular geometry from the impact 
energy E and the loss energy Tilo by 
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All experiments have been carried out in an UHV-system 
equipped with a HREELS spectrometer "Delta 0.5" orig- 
inally designed by Ibach et al. [13| ■ The pressure was kept 
below 1 x 10~ 10 Torr during all measurements. 

The graphene/silicon carbide sample used for this in- 
vestigation was prepared ex-situ by a hydrogen etching 
step and the following graphitisation under atmospheric 
argon pressure This resulted in a 1.5 ML graphene in 
addition to the buffer layer measured by XPS (not shown 
here)[5j. After transfer to the ultrahigh vacuum (UHV) 
system the sample showed almost no contamination with 
hydrocarbons or water. In contrast, the hydrogen etched 
silicon carbide samples used for comparison, shows small 
amounts of hydrocarbons and dissociated and/or non dis- 
sociated water [15| . 

Fig. [Tj represents an energy-loss spectrum in specu- 
lar direction for 6H-SiC(0001) before (a) and after (b) 
graphene formation. As already shown [l5| the sur- 
face optical phonons, also called Fuchs-Kliewer phonons 
[IS], are totally quenched. A two peak structure with 
maxima at lo + = 1270 cmT 1 and w_ = 560 cmT 1 re- 
sults, which is shifted to lower (higher) energy values for 
higher (lower) primary beam energies and correspond- 
ingly smaller (higher) -values according to equation [1] 
and our measured data, as shown in Fig. [2] (red dots). 

For the explanation of the dispersion shown in Fig. 
[2] (blue line) we consider a thin conducting layer with 
a charge carrier plasmon, brought on top of a sub- 
strate with strong surface optical phonons with loso = 
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These two modes couple, 
which can be described in the framework of dielectric 
theory [H, 0, [l5|, [l7| ■ The dielectric function of a phonon 
is given by 
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and the dielectric function of a volume charge carrier 
plasmon by 
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FIG. 2: Dispersion of the coupled phonon plasmon modes u±. 
The red dots represent our measurements, taken at different 
primary beam energies (2.5 eV < Eq < 80 eV), the blue 
lines result from theory (equation [5]) . The dashed line indi- 
cates a yjq behavior, which would be typical for a 2DEG. The 
gap between loso and loto is indicated by full horizonal lines 
(black). The shaded area represents single particle intraband 
excitations (SPEintra). 



where the plasmon frequency is a constant, which is de- 
termined by the charge carrier density and material prop- 
erties as discussed below. 

The description as a volume plasmon is not sufficent here, 
a two dimensional plasmon mode depending on the thick- 
ness d and showing a dispersion of cup2d(q\\) oc has 
to be taken into account: 
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Neglegting phonon and plasmon damping (7 = and 
r = 0) and restricting the calculation to the region where 
qud << 1 allows the determination of an exact expression 
for the dispersion of the coupled modes [|[ : 
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These coupled modes have mixed phonon and plasmon 
characteristics. Referring to Fig. [21 at low momentum 
transfer (gii — > 0) the oj+ mode behaves more like a 
phonon mode and converges to the frequency of the clas- 
sical surface optical phonon loso ■ The ui _ mode vanishes 
to zero in the long wavelength limit and behaves similar 
to a classical two dimensional plasmon mode. On the 
other hand, at high momentum transfers the w_ mode 
converges to the frequency of the TO phonon and the 
lo + mode behaves like a plasmon. It has to be noted, 
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FIG. 3: HREELS spectrum (open circles) taken at 20eV pri- 
mary beam energy together with simulation (red line). The 
inset gives details of the two layer model. T = 130if . 



that this anticrossing also exhibits an energetically for- 
bidden zone between the TO phonon frequency and the 
surface optical phonon frequency. 

Fig. [2] displays the resulting dispersion together with 
the theoretical dispersion of the two coupled modes, 
which follow from equation [5] A plasmon frequency of 
12000 cmT 1 was fitted and a TO phonon frequency of 
760 cm -1 was used [IJj]. Except for the w+ mode at 
high momentum transfers the curves fit perfectly to the 
measured data. The reason for the discrepancy at high 
momentum transfers lies in the non sufficient satisfaction 
of the condition g||<i << 1. 

In contrast to material systems investigated earlier 
(e.g. silver on GaAs) Q, in the graphene on silicon car- 
bide heterosystem this dispersion relation can be directly 
verified (see Fig. [5]). In other systems so far, the plasmon 
damping is always notably higher than in the graphene 
silicon carbide system. Therefore the dispersion relation 
calculated did not fit the dispersion observed in these 
systems. In particular, here no u>so mode is visible, 
which arises as a third solution in eq. [5l when taking 
finite damping into account. A closer look on the plas- 
mon damping can be gained by simulation of the whole 
HREELS spectrum. 

It is well known from dipole scattering theory [l3l [l9j|. 
that the energy loss probability in a HREELS experiment 
is given by 
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The integration limits of (?|| are determined by the angu- 
lar dimension of the spectrometer aperture. Convoluting 
this loss probability with a suitable spectrometer func- 
tion, also taking into account temperature effects and 



double losses allows us to directly compare the simula- 
tions with the measured spectra[20]. 

The effective dielectric function to describe the 
graphene silicon carbide heterosystem is build out of the 
elementary dielectric functions of the phonon and plas- 
mon contributions (see eqs. [5] and [3]). The model in- 
cludes two layers, one for the silicon carbide substrate 
and one for the graphene overlayer. The dielectric func- 
tion for the graphene overlayer just includes the charge 
carrier plasmon and is therefore identical to eq. [3] For 
the silicon carbide layer both, a charge carrier plas- 
mon from the slightly n-doped silicon carbide beneath 
(n 1.2 x 10 15 cm~ 3 ) and the TO phonon, have been 
included: 



2 

TO 



UJ 



PL 



UJ 



TO 



IJUJ 



ir 



SiCU 



(7) 



The effective dielectric function, which in contrast to the 
single layer dielectric functions now also takes the par- 
allel momentum transfer q» into account, can now be 
calculated [2lj|22]: 
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Although the basic theory behind this calculation is sim- 
ilar to the simple dispersion calculation shown above, it 
allows us to take the plasmon and phonon damping into 
account. 

Fig. [3] shows the numerical simulation (red line) to- 
gether with the corresponding HREELS-spectrum (open 
circles) . The simulation fits the measured data quite well, 
even though no band bending or interface effects (i.e. 
change of the charge carrier concentration in the silicon 
carbide as a function of depth) have been taken into ac- 
count. For other primary beam energies apart from the 
Eo = 20el/-spectrum shown in Fig. [3] simulations were 
carried out, too (not shown here). These simulations do 
fit the experimental results also very well. In general, 
the higher the primary beam energy is chosen, the better 
the simulation fits the experimental data. In addition 
it should be mentioned that with higher primary beam 
energies the uj + mode exhibits a Fano type resonance 
shape 23] . Further discussion of this behavior is beyond 
the scope of this contribution. 



TABLE I: Parameters for the simulation of the graphene sil- 
icon carbide heterostructure HREELS spectrum 
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Table U summarizes the parameters used in the simu- 
lation. The plasmon frequency is about 500 cmT 1 lower 
than in the dispersion simulation. This can be attributed 
to taking the damping into account. 

From angle resolved photoelectron spectroscopy it is 
well known that graphene in the graphene silicon car- 
bide heterosystem is electron doped. The measurements 
[j^ - Eij l show, that the linearly dispersing 7r-bands which 
provide the free charge carriers for the plasmons ob- 
served start about Ep = A50meV below the Fermi level 
of the heterosystem. Using the relations Ep = fivpkp 
and kp = ^Airn/ (g s g v ) (with t? , ? = 2 and g v = 2 for 
the spin and valley degeneracies) [27] , the sheet electron 
density is determined to be 1.5 x 10 13 cto~ 2 . If we con- 
sider a typical plasmon which is found in a two dimen- 
sional electron gas by setting wp/, = yj implying a 
parabolic component in the band structure, we obtain 
an effective mass of 0.03 x m e . The plasmon damping 
Tpp = e/(/zm*) can be used to calculate the electron 
mobility /i = 1700cm 2 /(Fs). These values fit the pic- 
ture of previously published values from Hall mobility 
measurements @. 

On the other hand, Das Sarma et al. [1(| calculated 
the plasmon frequency for graphene depending on the 
electron density for a two dimensional dirac system to 
be: 



UP2D = \ —z (47r-n 2 _ D ) / VFy/qT (9) 
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This model for the plasmon frequency takes the linear 
band dispersion in graphene into account. Using the ex- 
perimental data fit derived with eq. U] one can calculate 
the electron density to be 1.4 x 10 11 cm -2 . For this den- 
sity the Fermi energy shift would only be Ep — 45meV", 
one order of magnitude lower than measured. The doping 
of graphene on silicon carbide seems to be too high, so 
that a nonlinear contribution to the electronic dispersion 
has already to be considered. 

We have shown, that the strong coupling between the 
charge carrier plasmon in graphene and the surface op- 
tical phonon in silicon carbide can be understood in 
terms of dielectric theory. The dispersion differs strongly 
from the phonon and plasmon dispersion in an uncoupled 
case. This understanding of the coupled modes is espe- 
cially important for transport calculations in the epitax- 
ial graphene silicon carbide heterosystem. The scatter- 
ing processes with the coupled plasmon phonon modes, 
in particular in the low momentum low energy regime 
strongly depend on the changed dispersion of the w_ 
branch of the dispersion. 
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